The zero-dimensional O(N) vector model as a benchmark for perturbation theory, the 
large-iV expansion and the functional renormalisation group 
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We consider the zero-dimensional O(N) vector model as a simple example to calculate n-point 
correlation functions using perturbation theory, the large- iV expansion, and the functional renormal- 
isation group (FRG) . Comparing our findings with exact results, we show that perturbation theory 
breaks down for moderate interactions for all N, as one should expect. While the interaction- 
induced shift of the free energy and the self-energy are well described by the large- TV expansion even 
for small N, this is not the case for higher-order correlation functions. However, using the FRG in 
its one-particle irreducible formalism, we see that very few running couplings suffice to get accurate 
results for arbitrary N in the strong coupling regime, outperforming the large- iV expansion for small 
N. We further remark on how the derivative expansion, a well-known approximation strategy for 
the FRG, reduces to an exact method for the zero-dimensional O(N) vector model. 

PACS numbers: ll.10.Hi, ll.10.Jj, 11.15.Pg 



I. INTRODUCTION 

Field theory and its applications to condensed mat- 
ter and high-energy physics are difficult subjects in the 
natural sciences. The solution of almost any field the- 
ory requires approximations which rely on the mastery 
of complicated expressions such that it is easy to get 
bogged down in technical details. As there is no exact so- 
lution to compare with, the quality of results from differ- 
ent approximation strategies like perturbative methods, 
the large- N expansion, or non-perturbative methods such 
as the functional renormalisation group (FRG) is often 
difficult to judge. 

An instructive problem to learn about the conceptual 
underpinnings of asymptotic series (see ej*. Refs. [l|-y]) 
and Feynman diagrams (see e.g. Ref. |4J) in a simple 
context is the evaluation of the integral (with r real and 
u > 0) 



/°° / r u \ 

d^exp^--^ 2 --y + J^ 



(1.1) 



For u — Q and r > 0, this integral reduces to a Gaus- 
sian integral which can be calculated exactly, resulting 
in Z (J) = yjlitjr exp (j 2 /2r). Keeping u finite, Z(J) 
corresponds to the classical partition function of a par- 
ticle in an anharmonic potential. Alternatively, we can 
also think about this partition function as a toy model 
of a (zero-dimensional) field theory. By taking deriva- 
tives of Z{J) with respect to J, it is possible to obtain 
all n-point functions (also called correlation functions or 
Green functions), 



G (r 



1 d n Z{J) 



Z(0) dJ n 



(1.2) 
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Since the action S(<p) = \ip 2 
the Zi transformation ip — > 



h ttV i s invariant under 
-if, all n-point functions 



with odd n vanish. In fact, for the zero-dimensional field 
theory considered here, all n-point functions can also be 
calculated by taking successive derivatives of Z — Z(0) 
with respect to r, but as this feature does not generalise 
to a generic field theory, we follow here the usual strategy 
of including the source term Jtp and taking derivatives 
with respect to the sources. It turns out that the in- 
tegrals determining the n-point functions can all 
be evaluated analytically (see e.g. Ref. [l[ for the case 
n = 0). However, in order to learn about perturba- 
tive many-body techniques, it is instructive to expand 
Eq. (|1.1|) or Eq. (|1.2|) for r > in powers of u. As one 
should expect, such a perturbative approach is justified 
for small u/r 2 and breaks down for u/r 2 — 0(1). 

In order to go beyond simple perturbation theory one 
usually has to sum up an infinite series of diagrams. But 
which classes of diagrams should be included? What if 
there is no small parameter to expand in? A useful strat- 
egy is to duplicate an internal field (in our case (p) to an 
iV-component field and to calculate quantities in powers 
of \/N. This approach is called the large- N expansion. 
By including the first one or two terms, one hopes to ob- 
tain reliable results even for small N. Generalising the in- 
tegral in Eq. (jl.ljl to an A^-dimensional integral leads to a 
zero-dimensional A^-component field theory, which is also 
known as the zero-dimensional O(N) vector model 
As we will show below, the extrapolation of the leading- 
order large-A result to the cases N = 1 or vV = 2 gives 
rise to sizeable deviations from the exact result. While 
it is possible to derive higher-order coefficients in a series 
expansion in 1/N for the zero-dimensional O(N) vector 
model, this is usually problematic in more involved field 
theories. Clearly, there is a need for alternative meth- 
ods which allow for a non-perturbative solution of field 
theories. 

A non-perturbative method to determine all ?i-point 
functions is the functional renormalisation group (FRG). 
There are basically two different implementations of the 
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FRG (also called non-perturbative RG): The Wilson- 
Polchinski approach |6| and the effective average action 
approach [7-9]. While the Wilson- Polchinski approach 
was used with much success in proofs of renormalisabil- 
ity [1, [To| , the effective average action approach is better 
suited for most practical purposes [lll - ll9| . The reason 
for the success of the latter is grounded on the fact that 
it is based on the generating functional of one-particle 
irreducible vertices (also known as the effective action). 

The evaluation of the partition function Z (J) given in 
Eq. (jl.lj) has already served as a simple example for an 
application of the FRG [H M&k ( see also Ref - (H)- 
Within the FRG, a regulator is added to the Gaus- 
sian part of the action. This regulator is usually im- 
plemented in such a way that it gives an extra mass to 
low-energy fluctuations. For the zero-dimensional O(N) 
vector model it effectively replaces the mass term r by 
r + i?A, where initially R\ — > oo such that the dimension- 
less parameter u/(r + R\) 2 — > and the initial condition 
can be determined exactly. 

The central quantity to compute within the framework 
of the one-particle irreducible FRG is the effective action. 
This effective action is the generating functional of all 
one-particle irreducible vertices (for the zero-dimensional 
O(N) vector model all functionals reduce to functions). 
All n-point functions can be derived from this generating 
functional. In the presence of a regulator, the effective 
action is replaced by an effective average action which 
for i?A — > oo turns out to be equal to the action itself. 
Gradually removing the regulator, the flow of the effec- 
tive average action is governed by a flow equation 0-Q 
which is simple, but nevertheless exact. During the flow, 
an infinite number of couplings is generated, and trying 
to determine the exact solution is usually a hopeless en- 
deavour [2~ij . Finding a good truncation by keeping some 
couplings, but neglecting others, is maybe the most dif- 
ficult step in applying the FRG. However, as we will see, 
for the zero-dimensional O(N) vector model, already very 
few coupling constants can be sufficient to obtain a very 
good description of the strong-coupling limit. 

The remainder of this paper is organised as follows: In 
Sec. |TT] we generalise Eq. (II. lj) to the zero-dimensional 
O(N) vector model, introduce its generating functions 
and discuss its symmetries. After presenting an exact so- 
lution for all n-point functions in Sec. lIIIl we demonstrate 
how to calculate the n-point functions perturbatively in 
Sec. |IV] Some knowledge about the linked-cluster theo- 
rem and Wick's theorem are helpful here. We then dis- 
cuss different versions of the large- A expansion in Sec.fVl 
Higher-order results for both perturbation theory and the 
large-A expansion are deferred to the appendix. A brief 
introduction to the FRG is given in Sec. IVII ^25| . where 
we also study the zero-dimensional O(N) vector model 
by applying two widely used approximation strategies: 
While only three coupling constants suffice to give quite 
reliable results for the free energy and the two-point func- 
tion within the vertex expansion, for the simple problem 
considered here, the FRG flow equations can in princi- 



ple be calculated by directly solving a partial differen- 
tial equation numerically. In fact, another commonly 
used approximation strategy, the derivative expansion 
(which for the zero-dimensional O(N) vector model has 
no derivative term), reduces to an exact method. Finally, 
we discuss and compare our results in the conclusion in 
Sec. lVIll and make a few further remarks on how the stud- 
ies of this paper generalise to more involved field theories. 



II. GENERATING FUNCTIONS AND 
SYMMETRIES OF THE ZERO-DIMENSIONAL 
O(N) VECTOR MODEL 

Generalising Eq. to an A-dimensional integral 

and defining J T>(p := J\iLi I^oo ^Pii the generating 
function of all Green functions reads 

Z{J) = J V0 exp (-S(0) + J ■ <pj . (2.1) 

Here, 

S(0) = S o (tf+S int (0) = r -0 2 + ^(0 2 ) 2 (2.2) 

is a generalisation of the action given in Eq. (jl.ip and 
Cp 2 = 5Zj =1 <p 2 is the squared norm of the vector tp. While 
only Re (u) has to be positive for the integrals to con- 
verge, perturbation theory also requires Re (r) > 0. Even 
though both r and u can be considered to be complex, we 
will be mainly interested in the case where both of these 
parameters are real. However, let us mention here that 
the exact and perturbative results to be derived in the 
following two sections are also valid for complex r and u 
(provided Re (u) > for reasons of convergence and also 
Re (r) > in the perturbative case). For real r ^ 0, it 
is in principle possible to rescale tp — > <p/\/\r\, thereby 
replacing r by ±1 and u by the dimensionless parameter 
u/r 2 . However, let us not do so. We thereby avoid hav- 
ing to distinguish the cases r < 0, r = 0, and r > (in 
which we will be mainly interested). Furthermore, 1/r 
can be interpreted as the free propagator, and we would 
rather not replace this free propagator by one. 

The generating function Z(J) stores the information 
about all n-point functions, albeit for many purposes not 
in its most convenient form. The generating function 
of all connected Green functions is obtained by simply 
taking the logarithm of Eq. (|2.ip (for a neat proof based 
on the replica trick, see e.g. Ref. 



(for : 
0), 



W(J) = In 



Z{J) 
Z 



(2.3) 



The normalisation with the partition function Zq = 



Z(J = 0) 



(2ir/r) N / 2 of the free system lends W(0) 



the interpretation of the negative interaction-induced 
shift of the free energy. Calling the expectation value 
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of the field (p in the presence of sources 



(2.4) 



we can perform a Legendre transformation to make the 
transition from the generating function of connected 
Green functions to the generating function of one-particle 
irreducible vertices [2] , 



r(0) = j(4>) 



(2.5) 



This is similar to using the Legendre transformation to 
switch from the Helmholtz free energy F to the Gibbs 
potential G, which often turns out to be the more useful 
quantity. 

As a consequence of the generalisation of the simple 
integral in Eq. (jl.ip to N fields, an n-point Green func- 
tion now carries n indices ii,...,i n , ranging from 1 to 
N. Denoting the average with respect to the full action 
by 



fvgs 



(2.6) 



where nk is the number of indices i which are equal to 
k. It should be noted that all Green functions or vertices 
with an odd number of indices of the same value vanish. 
From now on, we will only consider := G[ n \, as 

it contains all the relevant information. For instance, it 
directly follows from Eq. (|2.9[) that 

G[% 2 = G%/3 = G^/3. (2.10) 

Since all other generating functions share the O(N) sym- 
metry, they can also be expanded like Z (J) in Eq. (12.81) . 

The relations between the different generating func- 
tions imply that their expansion coefficients are related. 
Taking derivatives of the above equations yields the fol- 
lowing useful identities after some straightforward ma- 
nipulations: 

#)=-r (0) , (2.ii) 



G*°> = In 

G (2) = G (2) = (Y(2) j * ^ (212) 
G (4) = G (4) _ 3 ^ G (2)^ 2 = _ |^ G (2)^ p(4) _ (213) 

Similar identities can also be derived for higher-order cor- 
relation functions. 



the n-point functions can be defined as 



G 



(n) 



{Vix ■■■fin) 



1 d n Z{J) 



■ dJ L 



(2.7) 



It turns out that for the present model, any (non- 
vanishing) Green function of order n already determines 
all other Green functions of the same order. To prove 
this statement, let us note that, by construction, our ac- 
tion given in Eq. (12.21) is invariant under O(N) rotations 
and therefore (as there are no derivative terms in zero 
dimensions) depends only on (p 1 . Since the action does 
not favour any particular direction in the iV-dimensional 
space, the same also has to be true for all generating 
functionals. We therefore have Z(J) = Z(J 2 ), W(J) = 
W(J 2 ), and T($) = f((f 2 ). A more formal proof of the 
fact that any symmetry of the action is also a symmetry 
of the generating functionals can be found, for example, 
in Refs. [13 ■ It follows that 



Z(J) = £ G< 



2m) 



(J^P 

(2m)! 



E^ ( 



2m) 



T N J 2 



(2m)! 



JV 



- e G(2m) (2m) , e [i i2 ™, . ; iN ) n j . 

m=0 v ' ii + ...+ijv=m v ' fc=l 



(2.8) 

and, by comparison with a generic Taylor expansion of 
Z(J), one finds 



N 



G 



(n)! AA(n fc /2)! 



(2.9) 



III. EXACT SOLUTION 

In order to calculate Z = Z(0), we use hyperspherical 
coordinates to obtain 



Z = J V0e~ s{ ^ = Sl N R N -i , 



(3.1) 



where FId is the surface area of the D-dimensional unit 
sphere and 



Rk — 



dire ! 4i =2 4 3 4 u 



fc + l fc+3 



+ 1\ _ fk + l 1 3r 2 

— J lFl l— ; 2 ; ^7 



(3.2) 



where iF\{a\b;z) is the Kummer confluent hypergeo- 
metric function. To determine G*- 2 * 1 = (<p 2 ), it is again 
helpful to exploit symmetry arguments and use (ip 2 ) = 

E^i^ 2 ) = N (<Pi) t0 arrive at 



G™ = -(p 2 ) = 

' NR N -! 



(3.3) 



Making use of Eq. (|2.10l) , one easily calculates 

{{0 2 ) 2 )=N{p\)+N{N-l){vW 2 ) 

_ N(N + 2) (4) 
3 ' 



(3.4) 
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FIG. 1. (colour online) Plots of £ for AT = 1. Within the 
resolution of this plot, the FRG results for n t = 6 lie on top 
of the exact results. 




FIG. 2. (colour online) Plots of r (0) for N = 2. Within the 
resolution of this plot, the FRG results for n t — 6 lie on top 
of the exact results. 



FIG. 3. (colour online) Plots of E for TV = 2. Within the 
resolution of this plot, the FRG results for n t — 6 lie on top 
of the exact results. 




FIG. 4. (colour online) Plots of T (4) for N = 2. The FRG 
results for n t — 6 lie slightly above the exact results. 



IV. PERTURBATION THEORY 



which results in 



N(N + 2) 



m 2 ) = 



N(N + 2) R N . 



(3-5) 



In principle, it is now also possible to go on to higher 
orders. 

Having found a method to compute the full Green func- 
tions in closed form, we can use Eqs. (|2.11j) - (|2.13j) 
to calculate Gc™^ and r( n ). A plot of the self-energy 
£ = r< 2 ) - r for N = 1 as well as plots of T^°\ E, and 
r^ 4 ) for N = 2 are shown in Figs. [TH4] as a function of 
the dimensionless parameter u/r 2 . In these figures, we 
also show perturbative and non-perturbative results to 
be discussed in the following sections. 



As a next step, we will take a look at calculating Green 
functions within perturbation theory. The usual strategy 
is to expand the interaction term exp (— ^j(</? 2 ) 2 ) in pow- 
ers of u. Introducing, in addition to the average with 
respect to the full action defined in Eq. (12.61) . the aver- 
age with respect to the Gaussian part So of the action, 



JV0e 



So(ip) 



\ / U --7 1 

the connected n-point functions are given by 

/ / -\ \ conn 

G(») = ^"e - *"*^ 

" N ' 1 / u\ conn 

\ i /ii _ _\ re i 



(4.1) 



fe=0 



(4.2) 
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The upper index "conn" indicates that only connected 
expressions should be considered when evaluating the 
Gaussian averages using Wick's theorem (see Refs. 0- 
4]). In fact, all occurring unconnected parts are cancelled 
by terms originating from an expansion of the partition 
function in the denominator of Eq. (|2.6[) , a result known 
as the linked- cluster theorem. 

In a nutshell, Wick's theorem states that 
((fi 1 . . . (pi 2m )o equals all possible contractions of pairs of 
fields, e.g. (<Pi 1 <Pi 3 <Pi 3 <Pi 4 )o = (^1^2)0(^3^4)0 + 
(^1^3)0(^2^4)0 + (fix <Pu)o(<Pi 2 <Pi 3 )o wi th 
(<Pi 1 <Pi a )o = Si lt i 2 /r. A large number of terms is 
generated and it is useful to identify the occurring terms 
with Feynman diagrams. In order to calculate the (D(u k ) 

contribution to Gc , one can apply the following 
Feynman rules: 

1. Draw all topologically inequivalent connected dia- 
grams Di with k wiggly 4-point vertices ^r^y^r^,^ 
and 2k + m propagators, which are denoted by a 
straight line . 

2. Assign an index 1 < j < N (or colour) to each line, 
keeping in mind that all external legs and both lines 
connected to each side of a vertex must have the 
same colour. 

3. For each diagram, determine its combinatorial fac- 
tor Ci, which is the number of Wick contractions 
it corresponds to (see below). 

4. For each diagram, count the number of colours that 
can be chosen independently and denote it by n.^. 

5. Assign a factor of \jr to each propagator and a 
factor of — it/4! to each vertex. Multiply the result 
by CiN ni /k\. 

6. The 0(u k ) contribution to G (2m) is the sum of all 
these diagrams. 

Given our Feynman rules, all that is left to do is to ap- 
ply them to determine G^ , G^ , and Gc . It is helpful 
to organise the diagrams by their number of vertices and 
number of independent colours. As many students have 
difficulties in figuring out the right combinatorial factors, 
let us discuss here the three examples depicted in Fig. [5] 

(a) To determine C a , consider the following reasoning: 
There are four locations to place the bubble, four 
ways of connecting Line 1 to the second vertex and 
two more ways of connecting Line 2 to the second ver- 
tex. Having done that, the location of the remaining 
line is already determined, hence C a =4x4x2 = 2 5 . 
As can be seen from the diagram, n c = 2 and there- 
fore its total contribution is 

1 f-u\ 2 1 , 2 N 2 u 2 

2UJ ^ X2 XN = 36^ ^ 




(a) C a = 2 5 (b) C b = 2 




FIG. 5. Examples of Feynman diagrams. The straight and 
wiggly lines denote the free propagator 1/r and the interaction 
u, respectively. 



(b) There are eight places to connect the first external 
leg to and four places for the second leg. Since this 
already fixes the internal line connecting the sides the 
legs are attached to, one has only an additional factor 
of two from arranging the lines that form the bubble 
and thus C b = 8 x 4 x 2 = 2 6 . It follows that the 
diagram contributes 

- — — x2 6 x N = . (4.4) 

2 V 4! J r 5 18 r 5 v ; 



(c) One has four choices for attaching the two lines at V3 
to a vertex, four choices for the first line connecting 
V3 to another vertex and two more choices how to 
connect V3 to V±. Since there are exactly 4! ways to 
shuffle around the external legs, C c = 4x4x2x4! = 
2 5 x 4!. Applying the remaining Feynman rules yields 

1 f-u\ 2 1 - , 2u 2 , s 
- — r -^x2 5 x4! = — =• . 4.5 

2 V 4! J r 6 3 r 6 v ; 



Adding up all the contributions for the connected 
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FIG. 6. Contributions to £ up to 0(u 2 



Green functions leads to 



N 2 + 2N u N 3 + 5N 2 + 6N ■ 



24 



144 



O u 3 



) = l_jV ± 2, + ^ + 5iV + 6^ +0 
r 6 r a 18 r 5 v y 



G « 1) = _» + 5« + 16^ +0( „, ) 



(4.6) 

(4.7) 
(4.8) 



With the connected Green functions determined, we can 
apply Eqs. (|2.11[) — (|2.1 3[) to finally calculate the irre- 
ducible vertices. Alternatively, we obtain the irreducible 
vertices by simply ignoring the unconnected diagrams, 
stripping off the external legs, and finally multiplying 
the result by —1. For n = 2, it is customary to subtract 
the inverse free propagator r from T^ 2 \ leading to the 
self-energy S = — r, whose graphical representation 
is shown in Fig. [5] Up to second order in u, we obtain 



fn i N 2 + 2N u N 3 + 5N 2 + 6N u 2 , , N 
r(°) = — r-. * tt: ~J + O (u 3 ) , 



24 r 2 



144 



(4.9) 



x = r(V-r = ^^- N2 + m + 8 K + Q(u 3 ) 



6 r 



:-!() 



r (4) 



N ■ 



6 r x 



+ O (u 3 ) 



(4.10) 
(4.11) 



Comparing the perturbative results for T^°\ E, and 
(as depicted in Figs.[THI]), calculated within perturbation 
theory up to sixth order (see appendix), with the exact 
solution, we see that perturbation theory is only reliable 
for small u/r 2 . As can be seen clearly in the figures, the 
usefulness of perturbation theory is restricted to small 
values of u/r 2 and going to higher order in perturbation 



theory does not help. This owes to the fact that our series 
expansions are not convergent at all and are only asymp- 
totic series 0-0] • Roughly speaking, the vanishing radius 
of convergence is a natural consequence of the fact that 
all integrals above are actually divergent for Re (u) < 0. 
More rigorously, it is possible to deduce this vanishing 
radius of convergence from the series of expansion coeffi- 
cients. Even though the series are not convergent, trun- 
cating the series expansions at a finite given order can 
give rise to excellent results for very small values of u/r 2 . 
In fact, for any given value of u/r 2 < 3/(2(7V — 1)), there 
is an optimal order 



3rj 
u 



(N-l) + 



9r 4 6(N - l)r 2 



(4-12) 

beyond which an inclusion of higher-order terms gives 
rise to less accurate results. While the residual error can 
be reduced for intermediate coupling strengths by using 
Borel summation techniques [2j, there clearly is a need 
for going beyond perturbation theory. 



V. LARGE- iV LIMIT 

The fact that N is an intrinsically dimensionless pa- 
rameter independent of any physical scale makes it es- 
pecially suitable for taking limits [26]. One expects that 
in the large- N limit fluctuations are suppressed due to 
the infinite number of degrees of freedom to couple to 
and therefore some aspects of our theory might actually 
simplify. To obtain a sensible limit, one demands that 
all terms appearing in the action [see Eq. (|2.2[) ] be ex- 
tensive. Since 1/r is just the free propagator, it does 
not scale with N, thereby implying that (p 2 — O(N). 
In order for S(ip) to be an extensive quantity, one thus 
finds that u — 0(1/ N) and it is reasonable to make the 
replacement u — > u/N. 

In the limit of large N, the leading and the next to 
leading order terms of the partition function Z can be 
calculated within the saddle point approximation [27| . 
Using hyperspherical coordinates and making the vari- 
able substitution y = ip 2 /N, the partition function turns 
into 



z = n N N N / 2 



dy 



h 2y 

where the function f(y) is defined by 



f(y) = -^y + ^y 2 ~^Hv) ■ 



This function has its minimum (for r > 0) at 




(5.1) 



(5.2) 



(5.3) 
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which in the limit u —¥ is given by yo = 1/r. Expanding 
f(y) around this minimum up to second order in y — y 
and performing the Gaussian integral we obtain 



z = n N N N ' 2 



2tt 



1/2 



-Nf( V0 ) 



1 + 0{ N 



with 



and 



f(yo) = -j- + \- ^ ln (^) , 



/"(yo) = \ - f- 



(5.4) 
(5.5) 

(5.6) 



which for u — > reduce to /(yo) = (1 + lnr)/2 and 
f"(yo) = r 2 /2. We can now either use Stirling's formula 
to find an asymptotic expression for HnN n / 2 or notice 
that when calculating the interaction-induced shift of the 

free energy, 1^°) = — In (jg^J , this term drops out. It im- 
mediately follows that (see also Ref. Q for an alternative 
derivation) 



r (o) = N 



In (ryn) 

4 4 2 v ' 



+ \ In (2 - ry ) 



oil- 

1 N 



(5.7) 



Higher-order terms in 1 /N can be obtained by including 
fluctuation corrections to the saddle-point approximation 
(see appendix) and are generally much more difficult to 
compute. 

As for the two-point function G c 2 ^ = G^ 2 \ using the 
saddle point approximation again, we have 



G^ = 



Jo 



^ we" 



Nf(y) 



JO 2y c 



Nf(y) 



(5.8) 



i.e. yo is just the expectation value of y in the limit 
N — > oo. Terms of arbitrary order in 1/JV for any n-point 
function of the zero-dimensional O(N) vector model can 
be obtained by taking successive derivatives of I^ ) with 
respect to r (see appendix). 

To get a better understanding of which diagrams have 
actually been summed up, let us now study what the 
highest power of N in a connected n-point function is. 
Assume that we are computing the fcth order contribu- 
tion to this Green function. There are exactly k vertices, 
each of which has four points to connect a propagator 
to. Of those 4fc points, only 4fc — n remain unoccupied 
after attaching the external legs, which means that at 
most 2k — n/2 colours can still be chosen independently 
of the colours of the external legs. Additionally, we are 
only interested in connected diagrams and one needs at 
least k — 1 propagators to connect k vertices. Since every 
propagator connecting two vertices further reduces the 




(a) Bubble diagram (b) Pincer diagram 

FIG. 7. Building blocks of Feynman diagrams in the large- iV 
limit. 



number of independent colours by one, we are left with 
k+l—n/2 colours that can be chosen freely, resulting in a 
contribution ~ ]\[k+i-n/2^ Taking [ n ^ account the cou- 
pling constant u/N of the vertices, one finally finds that 
the highest order contribution of N is N 1 ~ n ^ 2 . Hence, 



G (n) = O 



1 



JSj-n/2-1 



(5.9) 



such that lim 



JV- 
.(0) 



Ct 



(n) 



(2) 

for n > 2 while G c remains 



finite and Gc = O(N), as we have seen above. In the 
limit N — > oo, this implies that the disconnected Green 
functions factorise, i.e. for i =/= j we have 

(^ 2 )H^)<^) and (tf)=3«tf» a . (5.10) 

As discussed in Ref. [28[, this is precisely the behaviour 
that one should expect, since the connected contribu- 
tions to the correlation function can be identified as 
interaction-induced fluctuations which average out for 
large N. 

(2) 

We would now like to calculate G c diagrammatically 
in the large-A^ limit, including terms of all orders in 
u. From the above consideration it is clear that only 
the diagrams with a maximum number of independent 
colours contribute, namely those that can be constructed 
from the two basic elements in Fig. [71 All other elements 
consist of pieces in which a propagator connects both 
ends of a vertex, thereby reducing the number of 
independent colours. 

To determine all relevant diagrams, it is helpful to in- 
troduce shorthand notation using pairs of parentheses. 
One can express the bubble [Fig. [7] (a)] through an empty 
pair of parentheses "()", where it is implied that two 
neighbouring bubbles "()()" are connected to each other. 
Furthermore, we denote by "(D)" a pincer [Fig. [7J (b)], 
enclosing some arbitrary sub-diagram D. Examples of 
the notation are depicted in Fig. [51 Fortunately, the com- 
binatorial factor of diagrams constructed from those two 
elements only depends on the number of vertices and not 
on the internal structure. One thus finds the combinato- 
rial factor of a diagram with k vertices to be 



C k = 4 k k\ 



(5.11) 



With Ct determined, all that is left to do is count in 
how many ways one can arrange k pairs of parentheses 
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(b) ()() 



FIG. 8. The first three diagrams and their short-hand nota- 
tion. 




FIG. 9. Self-consistent Hartree approximation of the full 
propagator. While the free propagator 1 jr and the interaction 
u = u/N are again denoted by a straight line and a wiggly 
line, the Hartree propagator G c is denoted by a double line. 



to give a valid expression. The answer to this well-known 
combinatorial problem is given by the Catalan numbers 



Ck. — 



(2k)l 



(fc + l)!fc! 



(5.12) 



Putting everything together, we find that the u k con- 
tribution to the connected 2-point Green function is 

Mr \ iir 7 ) Ck and that in the large- N limit we have 



oo 



(2fc)! 



r S *!(* + !)! 



u 

6r 2 



which is in fact convergent for |u| < 3r 2 /2 and reproduces 
Eq. (53). 

Examining the above contributions, it turns out that 
the large-TV expansion of the TV-component field the- 
ory is identical to a self- consistent Hartree approxima- 
tion. Recalling that the Bubble diagram above is also 
known as the Hartree diagram, one can write down a 
self-consistency equation for G c , 



(5.14) 



whose graphical counterpart is displayed in Fig. |5] 
Solving the above quadratic equation, we obtain 




yo 



(5.15) 
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FIG. 10. (colour online) Plot of the self-energy E at fixed 
u/r 2 — 3 as a function of the continuous parameter N. Even 
though our model system has direct physical meaning only 
for integer TV (labelled by circles, triangles, and boxes), an 
evaluation of E is possible in all cases also for non-integer 
values of TV. 



which is the exact result for G^ 2 ' = Gc 2 "* in the limit 
TV — > oo and reproduces our previous results. By keep- 
ing all terms of order 0(1/N) we also include the Fock 
term and many other diagrams. As a consequence, G< 2 ) 
then agrees for small u/r 2 with perturbation theory (and 
the exact result), but, in contrast to perturbation theory, 
it also gives accurate results even for small TV and larger 
values of u/r 2 . This can be seen in Figs. Q] and [3] In 
Fig. [10] we further plot E for a fixed value of u/r 2 as a 
function of N, which we can assume to be a continuous 
parameter. While E and r^ ) are quite well described by 
the large-TV expansion (if we include sufficiently many 
terms, see Pigs.[TK5}, the results for the vertex T^ 4 ) show 
considerable deviations from the exact result (see Fig. HJ) - 
Just as perturbation theory is an asymptotic expansion in 
(5.13) the dimensionless parameter u/r 2 , the large- N expansion 
is an asymptotic expansion in the dimensionless param- 
eter 1/N. Although this asymptotic expansion is well 
behaved for both and E up to N — 1, the asymp- 
totic expansion for clearly breaks down for N w 2 
(see Fig. [UJ. 



VI. FRG STUDY 

Finally, let us consider the effective average action 
FRG approach to our toy model. An FRG study of the 
N = 1 zero-dimensional field theory was proposed by 
Schonhammer [20l ] and can be found in the lecture notes 
by Meden [H[ or Pawlowski HH (see also Ref. [H| for 
a study of the case N = 2). The same problem is dealt 
with in Exercise 7.1 of Ref. [17], in which the reader is 
furthermore invited to show that due to its one-particle 
irreducibility the effective average action approach gives 
rise to much better results than the Wilson-Polchinski 
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FIG. 11. (colour online) Plot of the vertex f (4) = AT (4) 
at fixed u/r 2 = Nu/r 2 = 3 as a function of the continuous 
parameter N. 



approach. 

As has been mentioned in the introduction, in the effec- 
tive average action approach a large artificial mass term 
(regulator) is initially added to all low-energy fluctua- 
tions. Replacing the action S by S + AS a, a useful reg- 
ulator term for the zero-dimensional O(N) vector model 
is 



As A ($o = £y 



(6.1) 



This term effectively replaces the parameter r in the ac- 
tion by r + i?A- Subtracting a similar mass term from the 
flowing generating functional of one-particle irreducible 
vertices, 

T A ($) = Ja($)-$-W (Ja(^)) - AS A ($) , (6.2) 
this effective average action satisfies the initial condition 
lim T A ($)=S($) . (6.3) 

Ultimately, we are interested in the problem without the 
artificial regulator term. It is the aim of the FRG to 
follow the flow of various couplings as this regulator term 
is gradually removed. For Ra ->0we then obtain 



lim r A (0) = r($ 

Ha— >u 



(6.4) 



As shown by Wetterich and others 7 9], the effective 
average action obeys the following flow equation, 



5Ar A (0) = -Tr 



where Ra is just Ra times the unit matrix and 
( d d 



ij 



(6.5) 
(6.6) 



If we are also interested in the interaction-induced shift 
of the free energy, T^°\ we have to subtract from the 
r.h.s. of Eq. (|6.5[) the same term with u set equal to zero. 

One way of solving Eq. (|6.5[) is to expand both sides of 
this equation with respect to the fields and compare their 
coefficients to derive flow equations for the expansion co- 
efficients of Ta (</>)• This is the so-called vertex expansion. 
Proceeding along this way, one obtains a system of in- 
finitely many coupled ordinary differential equations, of 
which we list the first four, 



dATf 



N 



a A rf = - 



+ 



BaTT 



(OaRa) (Ga - G ,a) , 

r A 4) (OaRa) g\ , 
rf (8 a Ra) g\ 

{Yf)\dARK) 

Tf (OaRa) G\ 



N + 


2 


6 




N + 


4 


10 




AM 


-8 


3 




N + 


6 


14 



a 3 



(6.7) 
(6.8) 

(6.9) 



+ (A + i4) rf rf (d A RA)G 3 A 



5A+ 130 



(2) 

Here, Ga = 1/(T A ' + Ra) is the flowing propagator and 
Go, a = !/(»' + Ra) is the corresponding non-interacting 
propagator. Graphical representations of these flow 
equations are depicted in Fig. Q21 

We truncate our vertex expansion by setting all vertices 
higher than a given order n t equal to their initial values. 
Choosing n t <E {2,4,6}, the resulting system of ordinary 
differential equations consists of very few equations and 
can easily be solved numerically. As the couplings 
with n > 2 are completely decoupled from T^ ' , we do not 
need to keep track of the latter if we are not interested in 
the flow of the interaction-induced shift of the free energy. 
Even when truncating the system of flow equations at 
n t = 6 we therefore only need to consider three flowing 
couplings. According to Eq. (|6.3I) . our initial values for 
the rl™' 1 are 



r A 4) ) 3 (5 A i? A )Gi . (6.10) 



rf = r 

Ao 



for n > 6 



3.11) 



To numerically integrate our set of flow equations, we use 
the following regulator, 



Ra = A" 



i/(rf 



Go 



A 



(6.12) 
A, and Ga 



which leads to OaRa — -T/A 2 , 

A -1 — r). For N = 1, the resulting equations 

i. da eh. Since Ra diverges for A 



agree with Refs. [T7J, |21| 
and vanishes for A — > r _1 , our flow 
to be integrated from A = to A = r~\ As can 
seen in Figs. UHll the values obtained by setting nt = 6 
already turn out to be very close to the exact result and 



Since Ra diverges for 
i ..„,. n..™ equations have 
1 As can be 
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FIG. 12. (colour online) Graphical representation of the flow 
equations for the first four running couplings. While the n- 
point vertices T A n ^ are denoted by filled circles with n external 
legs and the derivative with respect to A is denoted by a dot, 




the flowing propagator = ( F)^ 1 + Ra ) is represented 



by a thick line. The cross appearing exactly once in every 
graph on the r.h.s. simply denotes the derivative of the regu- 
lator Ra with respect to A. 



thereby reduce the possible benefit that could be gained 
by including higher orders. 

So far, we have restricted our discussion of the zero- 
dimensional O(N) vector model to the case of positive 
r. While the perturbative approach is only meaningful 
for such values of r, the large- iV expansion and the FRG 
study can readily be generalised to negative r. For the 
large- N expansion we note that for general r the mini- 
mum of f(y) is located at yo = 9£ ^s 



^gn(ryi+^-r 
which for negative r leads to slight modifications of all 
equations. As concerns the (FRG) vertex expansion, the 
extension to negative r is possible by expanding the effec- 
tive average action T A ((f) around its flowing minimum. 

Besides the vertex expansion, the derivative expansion 
is a powerful approximation for the FRG. Within this 
approximation, terms of arbitrary order are kept for a 
homogeneous field configuration. As there are no spatial 
variations in a zero-dimensional field theory, the deriva- 



II 2 

< 



-Ra = 0.5|r| 
H A = \r 
Ra = 2|r 
Ra = 5|r 
Ra = 10|r| 




FIG. 13. (colour online) Evolution of AUa(p = <t> 2 /2) for 
u/r 2 — 1 and r > [upper (red) curves] as well as u/r 2 — 1 
and r < [lower (blue) curves] . It should be observed that the 
final (physical) curve with _Ra = always has a single-valued 
minimum located at p — 0, owing to the fact that there is no 
spontaneous symmetry breaking in a system of finite (or even 
zero) volume. 



tive approximation actually becomes an exact method to 
solve the zero-dimensional O(N) vector model. Because 
Ta^) is only a function of \<j>\, it is customary to de- 
fine p = <f 2 /2 = ^2iLi4>i/2 and introduce the flowing 
effective potential as U\(p) = T\ ((f)). Using 



d 2 U A (p) 



3.13) 



and letting 4> point along one of the coordinate axes, we 
have 



d d 

[5 A Ra] ( -=®-=tf A (p) + R A 
o<p 0(f> 



N 



d A R As ^fd 2 U A 



d A R A 



1 



N - 1 



2 \U' A { P ) + 2 P U'i{p) + R A U' A (p) + R A , 

(6.14) 

Defining U A as the expansion coefficients of U A (p) as in 



u A ( P ) = Yl 



A p n 



(6.15) 



71=0 



one can compare the series of U A (p) and r A (|0|), 



N p (2n) 

1 A 



r A(H) = E(t)i(2pr = E 



u 



(n) 



A p n 



(6.16) 



n=0 



11 



to find that 

r^ 2n) = (2n-l)!!C/| n) . (6.17) 

Solving the partial differential equation (|6.14j) for Ua(p) 
by discretising p, we obtain results which are very close 
to the exact result. We can then use Eq. (|6.17[) to obtain 
the vertex functions that we were initially interested in. 
In Fig. [13] we show the evolution of the effective poten- 
tial Ua{p) for both a positive and a negative value of r. 
This evolution is especially interesting for negative r. At 
the beginning of the FRG flow, all fluctuations are ex- 
cluded, such that the minimum of the effective average 
potential is located at a finite value of p. Neglecting fluc- 
tuations would therefore lead (as in mean-field theory) to 
a symmetry-broken state by picking one vector cf> which 
satisfies p — <f) 2 /2. But we do have to include fluctua- 
tions, and this is what the FRG does. As all fluctuations 
are integrated out, the minimum of U\(p) is moved to 
p = 0, giving rise to a single, non-degenerate ground 
state with 4> — 0. This result is a natural consequence of 
the fact that there is no spontaneous symmetry break- 
ing in any system of finite (including zero) size. Let us 
also note that the effective potential becomes a convex 
function of 0, as one should expect for a function whose 
origin lies in a Legendre transformation. It is impossible 
to recover this feature within the vertex expansion. 

VII. CONCLUSIONS 

In order to conclude our comparison of different ap- 
proaches to the zero-dimensional O(N) vector model, let 
us come back to the plots in Figs.[TH4] The most obvious 
fact to observe is that perturbation theory consistently 
breaks down for u/r 2 — 0(1) (or, more accurately, for 
u/r 2 = 0(1) if TV > 1). Furthermore, it can be seen 
that including higher orders of the expansion parameter 
u only exacerbates the situation in the strong coupling 
regime — exactly the behaviour one would expect from an 
asymptotic series whose contributions have alternating 
signs. 

Moving on, we see that when including the first two 
subleading orders of the self-energy E, the large- TV expan- 
sion makes surprisingly good predictions for TV — 2 and 
even TV = 1. However, the large- TV expansion fails for the 
calculation of the four-point vertex T^ 4 ) for both TV = 1 
and TV = 2, and one needs to go to much larger values 
of TV (or very small values of u/r 2 ) to see an asymptotic 
convergence. As concerns more involved model systems, 
terms of order 0(l/iV 2 ) are usually very difficult to com- 
pute, and it is widely known that the large- TV expansion 
can fail, as it does for instance for half-odd integer spin 
antiferromagnets (29| . 



Last of all, let us come back to the FRG results. 
Within the vertex expansion, one truncates the system 
of flow equations at a given order and including more 
terms in this expansion consistently improves this result. 
As can be seen in Figs. [2] [4[ for the results with n t = 6 
(keeping all vertices up to the six-point vertex) it be- 
comes difficult to distinguish the FRG results from the 
exact values no matter how large the interaction strength 
u/r 2 is. As higher-order vertices are more directly af- 
fected by our truncations, the results decrease in their 
precision when going to higher-order vertices. Despite 
of that, one can clearly see that the FRG outperforms 
the other approaches with only three running couplings. 
For the zero-dimensional O(N) vector model considered 
here, the FRG flow equation can directly be solved us- 
ing numerical methods. This is not possible in higher 
dimensions, as the space of possible functions <p(x) be- 
comes just too big [30| | . However, keeping in a derivative 
expansion in addition to the effective potential U\(p) just 
the derivative term as already present in the action of a 
if A theory leads to a very useful ap pro ximation, the lo- 
cal potential approximation (LPA) jlH . HH, 53] ■ It then 
turns out, that while the effective action is still convex, 
its minimum is not single-valued. Instead, one obtains a 
flat plateau, giving rise to spontaneous symmetry break- 
ing. 

Summing up, wc find that perturbation theory only 
works for small values of u/r 2 . On the other hand, both 
the large- TV expansion and the functional renormalisation 
group can be used to explore non-perturbative aspects of 
our field theory. However, as the large- TV expansion fails 
to determine the four-point vertex for TV = 1 or 2, the 
FRG turns out to be most reliable for small TV. 



ACKNOWLEDGMENTS 

We thank Holger Gies, Aldo Isidori and Peter Kopietz 
for valuable discussions. This work was supported by the 
German Research Foundation DFG through the research 
group FOR 723. 

Appendix A: Higher-order perturbation theory and 
results for the large-TV expansion 

In this appendix we collect some results for perturba- 
tion theory up to sixth order and also give the next-to- 
next leading order results of the one-particle irreducible 
rt-point vertices in the large-TV expansion for n < 4. 
These expressions are used in Figs.HHHas well as Figs.[TUl 
and [TO 

Generalising the calculations leading to Eqs. (|4.9[) - 
(|4.11[) to higher order in u using Mathematica, we obtain 
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(0) _TV 2 + 2TV u TV 3 + 5TV 2 + 6TV /u\ 2 5TV 4 + 44TV 3 + 128TV 2 + 120TV / u \ 3 
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\r 2 



10368 Vr 2 
21TV 6 + 386TV 5 + 29007V 4 + 110007V 3 + 207127V 2 + 151207V 
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To obtain the next-to-next leading order term in the 
large- TV expansion of F, we can expand f(y) in the parti- 
tion function Eq. (|5.1[) around its minimum up to fourth 
order in y — y and also expand the prefactor l/y up to 
second order in y — y . Collecting all the relevant terms, 
Eq. (|5.4I) gets replaced by 



z =n N N N ^ 2 



2tt 



1/2 



Avlf"{yo) 

\2r 2 y 2 - 27ry 



-Nf(y ) 



16 



67V(2-ryo) 3 



(A4) 



Dividing Z by Zq and taking the logarithm we obtain for 
the interaction-induced free energy shift 



r (o) = N 



ryo 1 l , , s 
— - - A - - 2 Mry ) 

(8 + ry )(ry - l) 2 



1 



67V(2 - ry y 
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(A5) 



While higher-order terms in 1/TV can be evaluated by 
generalising the above approach, a simple identity be- 
tween various derivatives of the free energy allows for 
a recursive evaluation of higher-order terms Taking 
derivatives of r^ ^ with respect to r, we furthermore ob- 
tain 



1 

S = — 



2(1 -ryo) 



4(ry 



TVy (2 - 
l) 2 (3ry 



ryo) 
-1) 



N 2 y a (2 - ry ) 5 
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TV 3 2/ 2 (2-ry ) 7 



(A6) 
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